Código
# Librerías necesarias
require(ggplot2)
require(tidyverse)
require(ggfortify)
require(plotly)
require(kableExtra)
require(vroom)
require(leaflet)
library(lubridate)4 Visualización
Presentamos un ejemplo para visualización de datos de sismos de Argentina de los últimos 10 años.
Se llama sismo al proceso físico de liberación instantánea de energía acumulada en el interior de las rocas de la litósfera. La energía acumulada es producto de los esfuerzos tectónicos generados en las partes rígidas de la Tierra. Los sismos ocurren y han ocurrido históricamente, sin embargo, en la actualidad hay más registro de eventos sísmicos a causa de los avances tecnológicos que han permitido localizarlos y caracterizarlos con mayor exactitud.
El riesgo sísmico de una región es la probabilidad de consecuencias adversas ante un evento sísmico y resulta de la combinación entre la peligrosidad sísmica propia del lugar y de la vulnerabilidad de las construcciones.
El INPRES (Instituto Nacional para la Prevención Sísmica) afirma que no es posible, en el sentido amplio de la palarba, predecir un sismo (o terremoto, que en la jerga local alude a un sismo que ocasiona víctimas y daños severos). Se considera una predicción sísmica formal a aquella en la que se indica la fecha de ocurrencia, el lugar (coordenadas) de ocurrencia (con la profundidad) y el tamaño (magnitud), del evento por ocurrir. La sismología ha permitido conocer e identificar todas las regiones en el mundo que son potencialmente sísmicas, es decir, regiones donde siempre han ocurrido, ocurren y ocurrirán terremotos mientras la Tierra tenga vida. Aunque el tamaño de los terremotos que ocurrirán en alguna de estas regiones puede ser estimado de manera aproximada en función de datos históricos, sin embargo, la fecha del terremoto no puede ser determinada con exactitud.
Si bien hay estudios científicos sobre predicción sísmica que realizan países como China, Estados Unidos y Japón; la comunidad sismológica mundial en la actualidad, no ha logrado obtener, con rigor científico estricto, predicciones de terremotos.
A continuación, se muestran algunos posibles análisis para datos de sismos de Argentina durante los años 2012 a 2022, fuertemente relacionados con herramientas de visualización y descriptivas.
4.1 Exploración inicial
Los datos de sismos pueden obtenerse del buscador disponible en el INPRES, acá, en el que se registran eventos sísmicos en territorio argentino y diferentes variables asociadas a estos, en particular, si se trata de un sismo que se haya sentido o no.
Para lo que sigue, trabajamos con el dataset reducido sismos-arg.csv en el que se incluyen datos de sismos de diferente intensidad y magnitud en la región continental del país (excluyendo Tierra del Fuego) desde el 7 de enero de 2012, hasta el 18 de mayo de 2022.
El dataset sismos-arg.csv se encuentra acá y ha sido transformado a partir de datos scrappeados del buscador de sismos del INPRES. A continuación, se muestran 10 datos de dicho conjunto.
Código
knitr::kable(datos[1:10,], caption = "Dataset reducido")| Fecha | Hora | Latitud | Longitud | Provincia | Percibido | Magnitud | Profundidad |
|---|---|---|---|---|---|---|---|
| 2022-05-18 | 13:09:56 | -31.820 | -67.040 | San Juan | FALSE | 3.0 | 134 |
| 2022-05-18 | 10:18:34 | -24.216 | -67.139 | Salta | FALSE | 3.2 | 185 |
| 2022-05-18 | 09:40:09 | -31.262 | -68.708 | San Juan | FALSE | 3.3 | 104 |
| 2022-05-18 | 09:25:43 | -32.521 | -70.085 | San Juan | FALSE | 2.7 | 112 |
| 2022-05-18 | 08:45:59 | -28.192 | -67.376 | Catamarca | FALSE | 2.7 | 142 |
| 2022-05-18 | 08:42:33 | -32.043 | -70.030 | San Juan | FALSE | 2.7 | 114 |
| 2022-05-17 | 23:22:48 | -23.442 | -66.846 | Jujuy | FALSE | 3.6 | 210 |
| 2022-05-17 | 21:21:18 | -23.503 | -66.831 | Jujuy | FALSE | 3.6 | 210 |
| 2022-05-17 | 17:17:58 | -22.794 | -66.228 | Jujuy | FALSE | 4.0 | 273 |
| 2022-05-17 | 14:46:04 | -31.505 | -68.640 | San Juan | FALSE | 2.7 | 103 |
La variable ´Percibido´ es la que registra si la intensidad del sismo fue percibida, o no, por la población del territorio afectado. Los sismos son más frecuentes de lo que, quizás, supongamos. El siguiente mapa ofrece una posible visualización para los sismos percibidos y no percibidos de los últimos 10 años. Dentro del mapa, es posible activar la capa topográfica para relacionar, por ejemplo, los sismos con el relieve del territorio.
Código
sismos_arg <- datos
pal1 <-
colorNumeric(palette = c("deeppink2", "deeppink4"),
domain = sismos_arg$Profundidad)
pal2 <-
colorNumeric(palette = c("dodgerblue2", "dodgerblue4"),
domain = sismos_arg$Profundidad)
mapa_arg <- sismos_arg %>%
leaflet(options = leafletOptions(attributionControl = FALSE)) %>%
addTiles(group = "OSM (default)") %>%
addProviderTiles("OpenTopoMap",
group = "Topográfico",
options = providerTileOptions(opacity = 0.5)) %>%
addCircleMarkers(
data = sismos_arg %>%
filter(Percibido == FALSE),
lng = ~ Longitud,
lat = ~ Latitud,
fillOpacity = 0.15,
radius = ~ Magnitud,
color = ~ pal2(Profundidad),
stroke = FALSE,
group = "Sismos no percibidos"
) %>%
addCircleMarkers(
data = sismos_arg %>%
filter(Percibido == TRUE),
lng = ~ Longitud,
lat = ~ Latitud,
fillOpacity = 0.15,
radius = ~ Magnitud,
color = ~ pal1(Profundidad),
stroke = FALSE,
group = "Sismos percibidos"
) %>%
addLayersControl(
baseGroups = c("OSM (default)", "Topográfico"),
overlayGroups = c("Sismos no percibidos", "Sismos percibidos"),
options = layersControlOptions(collapsed = TRUE)
)
mapa_argEl radio de cada punto indica la magnitud del sismo, mientras que la intensidad del color, su profundidad. Los sismos también se clasifican de acuerdo con la profundidad de su foco: desde la corteza terreste y hasta 70 km de profundidad, son superficiales; entre los 70 y los 450 km de profundidad, son intermedios; y para más de 450 km profundidad, profundos. Los más profundos, en el mapa, se ubican en la capa de sismos no percibidos (azul) sobre la falla que recorre en línea más o menos vertical las provincias de Salta y Santiago del Estero. Profundidades intermedias se encuentran, mayoritariamente, en el oeste de Jujuy. El resto, en general, corresponde a profundidades superficiales.
Las provincias de San Juan, Salta, Jujuy, La Rioja, Mendoza, Catamarca, y Córdoba son aquellas que registran más de 1000 eventos durante el período observado. El siguiente mapa ofrece una visualización de sismos para dichas provincias, indicando con color su magnitud.
Código
top <-
c("San Juan",
"Salta",
"Jujuy",
"La Rioja",
"Mendoza",
"Catamarca",
"Córdoba")
pal3 <-
colorNumeric(
palette = c("gold", "sienna1", "firebrick", "darkred", "orangered4"),
domain = sismos_arg$Magnitud
)
mapa_arg_top <- sismos_arg %>%
leaflet(options = leafletOptions(attributionControl = FALSE)) %>%
addTiles(group = "OSM (default)") %>%
addProviderTiles("OpenTopoMap",
group = "Topográfico",
options = providerTileOptions(opacity = 0.5)) %>%
addCircles(
data = sismos_arg %>%
filter(Provincia == top[1]),
lng = ~ Longitud,
lat = ~ Latitud,
fillOpacity = 1,
color = ~ pal3(Magnitud),
stroke = FALSE,
group = top[1]
) %>%
addCircles(
data = sismos_arg %>%
filter(Provincia == top[2]),
lng = ~ Longitud,
lat = ~ Latitud,
fillOpacity = 1,
color = ~ pal3(Magnitud),
stroke = FALSE,
group = top[2]
) %>%
addCircles(
data = sismos_arg %>%
filter(Provincia == top[3]),
lng = ~ Longitud,
lat = ~ Latitud,
fillOpacity = 1,
color = ~ pal3(Magnitud),
stroke = FALSE,
group = top[3]
) %>%
addCircles(
data = sismos_arg %>%
filter(Provincia == top[4]),
lng = ~ Longitud,
lat = ~ Latitud,
fillOpacity = 1,
color = ~ pal3(Magnitud),
stroke = FALSE,
group = top[4]
) %>%
addCircles(
data = sismos_arg %>%
filter(Provincia == top[5]),
lng = ~ Longitud,
lat = ~ Latitud,
fillOpacity = 1,
color = ~ pal3(Magnitud),
stroke = FALSE,
group = top[5]
) %>%
addCircles(
data = sismos_arg %>%
filter(Provincia == top[6]),
lng = ~ Longitud,
lat = ~ Latitud,
fillOpacity = 1,
color = ~ pal3(Magnitud),
stroke = FALSE,
group = top[6]
) %>%
addCircles(
data = sismos_arg %>%
filter(Provincia == top[7]),
lng = ~ Longitud,
lat = ~ Latitud,
fillOpacity = 1,
color = ~ pal3(Magnitud),
stroke = FALSE,
group = top[7]
) %>%
addLayersControl(
baseGroups = c("OSM (default)", "Topográfico"),
overlayGroups = top,
options = layersControlOptions(collapsed = TRUE)
)
mapa_arg_topLa percepción del sismo por parte de la población no está únicamente asociada a la magnitud y/o profundidad. Por ejemplo, un sismo de magnitud 2.6 y profundidad de 7 km en Córdoba es percibido, mientras que otro, de magnitud 2.7 y profundidad de 12 km en Catamarca no es percibido. Lo mismo ocurre para magnitudes y/o profundidades mayores. El siguiente mapa ofrece una visualización similar al anterior, pero restringido únicamente a los sismos percibidos del conjunto de datos.
Código
mapa_arg_top_perc <- sismos_arg %>%
leaflet(options = leafletOptions(attributionControl = FALSE)) %>%
addTiles(group = "OSM (default)") %>%
addProviderTiles("OpenTopoMap",
group = "Topográfico",
options = providerTileOptions(opacity = 0.5)) %>%
addCircles(
data = sismos_arg %>%
filter(Provincia == top[1],
Percibido == TRUE),
lng = ~ Longitud,
lat = ~ Latitud,
fillOpacity = 1,
color = ~ pal3(Magnitud),
stroke = FALSE,
group = top[1]
) %>%
addCircles(
data = sismos_arg %>%
filter(Provincia == top[2],
Percibido == TRUE),
lng = ~ Longitud,
lat = ~ Latitud,
fillOpacity = 1,
color = ~ pal3(Magnitud),
stroke = FALSE,
group = top[2]
) %>%
addCircles(
data = sismos_arg %>%
filter(Provincia == top[3],
Percibido == TRUE),
lng = ~ Longitud,
lat = ~ Latitud,
fillOpacity = 1,
color = ~ pal3(Magnitud),
stroke = FALSE,
group = top[3]
) %>%
addCircles(
data = sismos_arg %>%
filter(Provincia == top[4],
Percibido == TRUE),
lng = ~ Longitud,
lat = ~ Latitud,
fillOpacity = 1,
color = ~ pal3(Magnitud),
stroke = FALSE,
group = top[4]
) %>%
addCircles(
data = sismos_arg %>%
filter(Provincia == top[5],
Percibido == TRUE),
lng = ~ Longitud,
lat = ~ Latitud,
fillOpacity = 1,
color = ~ pal3(Magnitud),
stroke = FALSE,
group = top[5]
) %>%
addCircles(
data = sismos_arg %>%
filter(Provincia == top[6],
Percibido == TRUE),
lng = ~ Longitud,
lat = ~ Latitud,
fillOpacity = 1,
color = ~ pal3(Magnitud),
stroke = FALSE,
group = top[6]
) %>%
addCircles(
data = sismos_arg %>%
filter(Provincia == top[7],
Percibido == TRUE),
lng = ~ Longitud,
lat = ~ Latitud,
fillOpacity = 1,
color = ~ pal3(Magnitud),
stroke = FALSE,
group = top[7]
) %>%
addLayersControl(
baseGroups = c("OSM (default)", "Topográfico"),
overlayGroups = top,
options = layersControlOptions(collapsed = TRUE)
)
mapa_arg_top_percAllí puede verse que, mientras que en San Juan los sismos percibidos parecen ser de magnitudes altas, en Córdoba ocurre lo contrario.
Por otro lado, la profundidad de los sismos guarda estrecha relación con la configuración geológica y tecnótica de su ubicación geográfica. El gráfico que sigue toma los eventos de las provincias con más registros y exhibe la relación entre la profundidad del sismo y los paralelos de su ubicación de registro. Llamativo resulta ver que la actividad sísmica disminuye en las profundidades que separan los sismos superficiales de los intermedios. Como veíamos, las mayores profundidades se detectan en las provincias de Salta y Jujuy, con la excepción de una pequeña región entre paralelos en San Juan.
Conocimiento sobre la configuración tectónica permite interpretaciones más profundas. Por ejemplo, la zona norte de Mendoza ubicada entre los paralelos -32.5 y -33.5 y donde se ubica parte del Gran Mendoza, se encuentra en la región de mayor peligrosidad sísmica según INPRES. Lo que ocurre es que allí se da la convergencia entre la placa Sudamericana y la placa de Nazca, que se extiende a lo largo de la fosa Peruano-Chilena y en su movimiento en la dirección Este, subduce debajo de la Sudamericana, que se mueve en sentido opuesto. Dicho mecanismo induce un régimen de esfuerzos compresionales a la placa Sudamericana y que resultan responsables de la orientación y sentido del desplazamiento de las fallas lo que, como consecuencia, caracteriza la sismicidad de la zona.
El siguiente gráfico muestra la profundidad en relación con la longitud de la ubicación del foco para la zona norte de Mendoza ubicada entre los paralelos -32.5° y -33.5°, zona de subducción entre la placa de Nazca y la Sudamericana (zona de Benioff). Un análisis experto puede halar en el gráfico evidencia de la placa de Nazca subducida y moviéndose en dirección Este respecto a la placa Sudamericana (más detalle sobre esto puede consultarse acá).
La tabla que sigue resume, por provincia, la cantidad de sismos registrados, los percibidos, así como la magnitud y profundidad mediana. De las provincias con mayor cantidad de eventos registrados, Jujuy y Salta registran los sismos más profundos, mientras que Córdoba registra sismos superficiales, con una intensidad mediana de 2.8, y es de las provincias con mayor tasa de sismos percibidos de entre las de mayor cantidad de eventos registrados. San Juan registra muchos sismos, más de la mitad de ellos de profundidad media a profunda, pero con una tasa de percepción relativamente baja, como Salta y Jujuy.
Código
knitr::kable(sismos_prov, caption = "Resumen por provincias")| Provincia | n | nperc | Magmed | Profmed |
|---|---|---|---|---|
| Buenos Aires | 7 | 5 | 3.7 | 31 |
| Catamarca | 3317 | 109 | 2.8 | 139 |
| Chaco | 3 | 0 | 3.3 | 231 |
| Chubut | 1 | 0 | 3.8 | 169 |
| Corrientes | 2 | 0 | 3.7 | 267 |
| Córdoba | 1145 | 313 | 2.8 | 21 |
| Entre Ríos | 1 | 0 | 3.8 | 32 |
| Formosa | 3 | 0 | 3.8 | 342 |
| Jujuy | 5270 | 64 | 3.2 | 223 |
| La Pampa | 5 | 4 | 3.6 | 22 |
| La Rioja | 4616 | 166 | 2.8 | 114 |
| Mendoza | 4319 | 303 | 2.7 | 31 |
| Neuquén | 231 | 66 | 3.0 | 8 |
| Río Negro | 1 | 0 | 2.9 | 98 |
| Salta | 5781 | 61 | 3.2 | 195 |
| San Juan | 29920 | 619 | 2.7 | 106 |
| San Luis | 568 | 69 | 2.8 | 22 |
| Santa Cruz | 46 | 22 | 3.0 | 14 |
| Santiago del Estero | 255 | 59 | 3.2 | 30 |
| Tucumán | 326 | 46 | 2.8 | 27 |
Código
knitr::kable(sismos_prov_porc, caption = "Porcentaje de sismos percibidos para casos de más de 1000 eventos")| Provincia | pperc |
|---|---|
| Catamarca | 3.286102 |
| Córdoba | 27.336244 |
| Jujuy | 1.214421 |
| La Rioja | 3.596187 |
| Mendoza | 7.015513 |
| Salta | 1.055181 |
| San Juan | 2.068850 |
4.2 Análisis para el caso de San Juan
San Juan registra casi 30.000 eventos para el período considerado de poco más de 10 años, con una mediana de 8 sismos diarios. Un histograma para la cantidad diaria de sismos se observa a continuación.
En el histograma puede verse que hay por lo menos una observación que supera los 100 sismos diarios. En efecto, se trata del terremoto del 18 de enero de 2021 en La Rinconada, San Juan, en el que se reportaron colapsos de viviendas de adobe, grietas distensivas a lo largo del eje de la Ruta Nacional 40, derrumbes parciales, y otros daños serios principalmente en los Departamentos de Pocito y Sarmiento. Durante las 72 hs posteriores al terremoto de magnitud 6.4, se registraron más de 150 réplicas de magnitudes entre 2.5 y 5.3. Más información acá.
El siguiente gráfico muestra la cantidad de sismos diarios. Puede apreciarse las réplicas de tal terremoto en las observaciones que corresponden a los días entre el 19 y el 23 de enero de 2021.
Como vimos, San Juan registra más de 600 sismos percibidos durante el período, algo más que el 2% del total de sismos registrados. La frecuencia con la que se perciben tales sismos es alta. Eso puede observarse en el histograma para la cantidad de días transcurridos entre un sismo percibido y otro.
Córdoba y Mendoza, registrando cerca del 27% y 7% de sismos percibidos respecto del total, respectivamente, tienen una frecuencia ligeramente más espaciada que la de San Juan.
4.2.1 Gutenberg-Richter
La ley de Gutenberg-Richter (1958) es una fórmula que permite cuantificar la relación entre la frecuencia y la magnitud de la actividad sísmica de una región. La relación fue inicialmente propuesta por Charles Francis Richter y Beno Gutenberg en 1954 y es aun objeto de caracterización de sismos, pues se ha mostrado relativamente invariante por región y en el tiempo.
N es un número de sismos de magnitud mayor que M durante un período de tiempo “t”, logN es un logaritmo del número de sismos con magnitud M La ley establece que la recurrencia sísmica de magnitud mayor o igual a \(M\) durante un período de tiempo \(t\), en una región, se relaciona con la cantidad de sismos, \(N\), de magnitud mayor que \(M\) durante el período, a través de
\[ \log _{10} N=a-b M, \]
donde \(a\) y \(b\) son constantes determinadas por la naturaleza sísmica de la región, en general, a partir de datos históricos.
Ajustamos los datos de San Juan al modelo de Gutenberg-Richter para estimar las constantes \(a\) y \(b\). En el siguiente gráfico puede verse la cantidad \(N\) de sismos de magnitud mayor que \(M\) por año (escala logarítmica en base 10). Como los años 2012 y 2022 no están completos para el año calendario, reducimos los datos al período 2013-2021.
La estimación por mínimos cuadrados de \(a\) y de \(b\) puede verse a continuación.
Los coeficientes estimados son \(6.1812\) para \(a\), y \(1.1433\) para \(b\). El coeficiente \(b\) suele ser aproximadamente igual a \(1\) para sismos tectónicos. Como los datos con los que se estimó son de frecuencia anual, la cantidad \(a/b\) puede interpretarse como, en promedio y una vez al año, la ocurrencia de un terremoto de magnitud \(a/b\) o superior. En este caso, \(a/b=5.4065\), lo que sugiere que, en promedio, una vez al año se registra en San Juan un terremoto de magnitud aproximada \(5.4\) o superior.
En el siguiente gráfico se superpone, en azul, la curva estimada de Gutenberg-Richter
\[ \log _{10} N=6.1812-1.1433 M, \]
a partir de los datos de San Juan.
4.2.2 Estimación de la probabiblidad de ocurrencia de futuros sismos
Del análisis de los resultados de la tasa o frecuencia media de ocurrencia por unidad de tiempo, suele derivarse la estimación de la probabilidad de ocurrencia de futuros eventos sísmicos en la región estudiada, bajo el supuesto fundamental de que el nivel de actividad sísmica de los últimos años se mantendrá aproximadamente igual, en promedio, para los próximos años.
Para evaluar dicho riesgo sísmico, se modela la probabilidad con una distribución de Poisson, pues se supone que los eventos sísmicos se producen de un modo aleatorio e independiente, donde los tiempos de origen, las coordenadas de los focos y las magnitudes son variables independientes entre sí. Es decir, que la ocurrencia de un evento sísmico no tiene influencia en la ocurrencia de otro, y que la probabilidad de que dos eventos sucedan en el mismo sitio y al mismo tiempo es casi nula.
En general, tales hipótesis no concuerdan completamente con la naturaleza física de la acumulación de energía en zonas sísmicas. Sin embargo, y a pesar de estas limitaciones, este modelo provee una medida del riesgo sísmico relativamente sencilla, interpretable y típicamente usada por la comunidad experta para proveer una cuantificación de riesgo aceptable para las regiones.
Entonces, bajo este modelo, la probabilidad de ocurrencia de al menos un sismo de magnitud mayor que \(M\) en los próximos \(t\) años se calcula como \[ P(t)=1-e^{-N_1(M) \times t} . \] donde \(N_1(M)=10^{(a-b M \log t)}\). Luego, el número promedio de años para que un sismo de magnitud mayor que \(M\) ocurra está dado por \[ T_R=\frac{1}{N_1(M)}. \] Tomando como referencia la estimación de \(a\) y \(b\) obtenida para San Juan, la estimación de la probabilidad de ocurrencia de al menos un sismo de magnitud mayor que \(M\) en San Juan en los próximos \(t\) años, para \(t\in\{1,5,10,20\}\), se observa en el siguiente gráfico.
Magnitud t_1 t_5 t_10 t_20
1 3.00 1.00000000 1.0000000 1.0000000 1.0000000
2 3.01 1.00000000 1.0000000 1.0000000 1.0000000
3 3.02 1.00000000 1.0000000 1.0000000 1.0000000
4 3.03 1.00000000 1.0000000 1.0000000 1.0000000
5 3.04 1.00000000 1.0000000 1.0000000 1.0000000
6 3.05 1.00000000 1.0000000 1.0000000 1.0000000
7 3.06 1.00000000 1.0000000 1.0000000 1.0000000
8 3.07 1.00000000 1.0000000 1.0000000 1.0000000
9 3.08 1.00000000 1.0000000 1.0000000 1.0000000
10 3.09 1.00000000 1.0000000 1.0000000 1.0000000
11 3.10 1.00000000 1.0000000 1.0000000 1.0000000
12 3.11 1.00000000 1.0000000 1.0000000 1.0000000
13 3.12 1.00000000 1.0000000 1.0000000 1.0000000
14 3.13 1.00000000 1.0000000 1.0000000 1.0000000
15 3.14 1.00000000 1.0000000 1.0000000 1.0000000
16 3.15 1.00000000 1.0000000 1.0000000 1.0000000
17 3.16 1.00000000 1.0000000 1.0000000 1.0000000
18 3.17 1.00000000 1.0000000 1.0000000 1.0000000
19 3.18 1.00000000 1.0000000 1.0000000 1.0000000
20 3.19 1.00000000 1.0000000 1.0000000 1.0000000
21 3.20 1.00000000 1.0000000 1.0000000 1.0000000
22 3.21 1.00000000 1.0000000 1.0000000 1.0000000
23 3.22 1.00000000 1.0000000 1.0000000 1.0000000
24 3.23 1.00000000 1.0000000 1.0000000 1.0000000
25 3.24 1.00000000 1.0000000 1.0000000 1.0000000
26 3.25 1.00000000 1.0000000 1.0000000 1.0000000
27 3.26 1.00000000 1.0000000 1.0000000 1.0000000
28 3.27 1.00000000 1.0000000 1.0000000 1.0000000
29 3.28 1.00000000 1.0000000 1.0000000 1.0000000
30 3.29 1.00000000 1.0000000 1.0000000 1.0000000
31 3.30 1.00000000 1.0000000 1.0000000 1.0000000
32 3.31 1.00000000 1.0000000 1.0000000 1.0000000
33 3.32 1.00000000 1.0000000 1.0000000 1.0000000
34 3.33 1.00000000 1.0000000 1.0000000 1.0000000
35 3.34 1.00000000 1.0000000 1.0000000 1.0000000
36 3.35 1.00000000 1.0000000 1.0000000 1.0000000
37 3.36 1.00000000 1.0000000 1.0000000 1.0000000
38 3.37 1.00000000 1.0000000 1.0000000 1.0000000
39 3.38 1.00000000 1.0000000 1.0000000 1.0000000
40 3.39 1.00000000 1.0000000 1.0000000 1.0000000
41 3.40 1.00000000 1.0000000 1.0000000 1.0000000
42 3.41 1.00000000 1.0000000 1.0000000 1.0000000
43 3.42 1.00000000 1.0000000 1.0000000 1.0000000
44 3.43 1.00000000 1.0000000 1.0000000 1.0000000
45 3.44 1.00000000 1.0000000 1.0000000 1.0000000
46 3.45 1.00000000 1.0000000 1.0000000 1.0000000
47 3.46 1.00000000 1.0000000 1.0000000 1.0000000
48 3.47 1.00000000 1.0000000 1.0000000 1.0000000
49 3.48 1.00000000 1.0000000 1.0000000 1.0000000
50 3.49 1.00000000 1.0000000 1.0000000 1.0000000
51 3.50 1.00000000 1.0000000 1.0000000 1.0000000
52 3.51 1.00000000 1.0000000 1.0000000 1.0000000
53 3.52 1.00000000 1.0000000 1.0000000 1.0000000
54 3.53 1.00000000 1.0000000 1.0000000 1.0000000
55 3.54 1.00000000 1.0000000 1.0000000 1.0000000
56 3.55 1.00000000 1.0000000 1.0000000 1.0000000
57 3.56 1.00000000 1.0000000 1.0000000 1.0000000
58 3.57 1.00000000 1.0000000 1.0000000 1.0000000
59 3.58 1.00000000 1.0000000 1.0000000 1.0000000
60 3.59 1.00000000 1.0000000 1.0000000 1.0000000
61 3.60 1.00000000 1.0000000 1.0000000 1.0000000
62 3.61 1.00000000 1.0000000 1.0000000 1.0000000
63 3.62 1.00000000 1.0000000 1.0000000 1.0000000
64 3.63 1.00000000 1.0000000 1.0000000 1.0000000
65 3.64 1.00000000 1.0000000 1.0000000 1.0000000
66 3.65 1.00000000 1.0000000 1.0000000 1.0000000
67 3.66 1.00000000 1.0000000 1.0000000 1.0000000
68 3.67 1.00000000 1.0000000 1.0000000 1.0000000
69 3.68 1.00000000 1.0000000 1.0000000 1.0000000
70 3.69 1.00000000 1.0000000 1.0000000 1.0000000
71 3.70 1.00000000 1.0000000 1.0000000 1.0000000
72 3.71 1.00000000 1.0000000 1.0000000 1.0000000
73 3.72 1.00000000 1.0000000 1.0000000 1.0000000
74 3.73 1.00000000 1.0000000 1.0000000 1.0000000
75 3.74 1.00000000 1.0000000 1.0000000 1.0000000
76 3.75 1.00000000 1.0000000 1.0000000 1.0000000
77 3.76 1.00000000 1.0000000 1.0000000 1.0000000
78 3.77 1.00000000 1.0000000 1.0000000 1.0000000
79 3.78 1.00000000 1.0000000 1.0000000 1.0000000
80 3.79 1.00000000 1.0000000 1.0000000 1.0000000
81 3.80 1.00000000 1.0000000 1.0000000 1.0000000
82 3.81 1.00000000 1.0000000 1.0000000 1.0000000
83 3.82 1.00000000 1.0000000 1.0000000 1.0000000
84 3.83 1.00000000 1.0000000 1.0000000 1.0000000
85 3.84 1.00000000 1.0000000 1.0000000 1.0000000
86 3.85 1.00000000 1.0000000 1.0000000 1.0000000
87 3.86 1.00000000 1.0000000 1.0000000 1.0000000
88 3.87 1.00000000 1.0000000 1.0000000 1.0000000
89 3.88 1.00000000 1.0000000 1.0000000 1.0000000
90 3.89 1.00000000 1.0000000 1.0000000 1.0000000
91 3.90 1.00000000 1.0000000 1.0000000 1.0000000
92 3.91 1.00000000 1.0000000 1.0000000 1.0000000
93 3.92 1.00000000 1.0000000 1.0000000 1.0000000
94 3.93 1.00000000 1.0000000 1.0000000 1.0000000
95 3.94 1.00000000 1.0000000 1.0000000 1.0000000
96 3.95 1.00000000 1.0000000 1.0000000 1.0000000
97 3.96 1.00000000 1.0000000 1.0000000 1.0000000
98 3.97 1.00000000 1.0000000 1.0000000 1.0000000
99 3.98 1.00000000 1.0000000 1.0000000 1.0000000
100 3.99 1.00000000 1.0000000 1.0000000 1.0000000
101 4.00 1.00000000 1.0000000 1.0000000 1.0000000
102 4.01 1.00000000 1.0000000 1.0000000 1.0000000
103 4.02 1.00000000 1.0000000 1.0000000 1.0000000
104 4.03 1.00000000 1.0000000 1.0000000 1.0000000
105 4.04 1.00000000 1.0000000 1.0000000 1.0000000
106 4.05 1.00000000 1.0000000 1.0000000 1.0000000
107 4.06 1.00000000 1.0000000 1.0000000 1.0000000
108 4.07 1.00000000 1.0000000 1.0000000 1.0000000
109 4.08 1.00000000 1.0000000 1.0000000 1.0000000
110 4.09 1.00000000 1.0000000 1.0000000 1.0000000
111 4.10 1.00000000 1.0000000 1.0000000 1.0000000
112 4.11 1.00000000 1.0000000 1.0000000 1.0000000
113 4.12 1.00000000 1.0000000 1.0000000 1.0000000
114 4.13 1.00000000 1.0000000 1.0000000 1.0000000
115 4.14 1.00000000 1.0000000 1.0000000 1.0000000
116 4.15 1.00000000 1.0000000 1.0000000 1.0000000
117 4.16 1.00000000 1.0000000 1.0000000 1.0000000
118 4.17 1.00000000 1.0000000 1.0000000 1.0000000
119 4.18 1.00000000 1.0000000 1.0000000 1.0000000
120 4.19 1.00000000 1.0000000 1.0000000 1.0000000
121 4.20 1.00000000 1.0000000 1.0000000 1.0000000
122 4.21 1.00000000 1.0000000 1.0000000 1.0000000
123 4.22 1.00000000 1.0000000 1.0000000 1.0000000
124 4.23 1.00000000 1.0000000 1.0000000 1.0000000
125 4.24 1.00000000 1.0000000 1.0000000 1.0000000
126 4.25 1.00000000 1.0000000 1.0000000 1.0000000
127 4.26 1.00000000 1.0000000 1.0000000 1.0000000
128 4.27 1.00000000 1.0000000 1.0000000 1.0000000
129 4.28 1.00000000 1.0000000 1.0000000 1.0000000
130 4.29 0.99999999 1.0000000 1.0000000 1.0000000
131 4.30 0.99999999 1.0000000 1.0000000 1.0000000
132 4.31 0.99999998 1.0000000 1.0000000 1.0000000
133 4.32 0.99999997 1.0000000 1.0000000 1.0000000
134 4.33 0.99999996 1.0000000 1.0000000 1.0000000
135 4.34 0.99999994 1.0000000 1.0000000 1.0000000
136 4.35 0.99999990 1.0000000 1.0000000 1.0000000
137 4.36 0.99999985 1.0000000 1.0000000 1.0000000
138 4.37 0.99999978 1.0000000 1.0000000 1.0000000
139 4.38 0.99999967 1.0000000 1.0000000 1.0000000
140 4.39 0.99999951 1.0000000 1.0000000 1.0000000
141 4.40 0.99999928 1.0000000 1.0000000 1.0000000
142 4.41 0.99999896 1.0000000 1.0000000 1.0000000
143 4.42 0.99999852 1.0000000 1.0000000 1.0000000
144 4.43 0.99999790 1.0000000 1.0000000 1.0000000
145 4.44 0.99999705 1.0000000 1.0000000 1.0000000
146 4.45 0.99999589 1.0000000 1.0000000 1.0000000
147 4.46 0.99999433 1.0000000 1.0000000 1.0000000
148 4.47 0.99999224 1.0000000 1.0000000 1.0000000
149 4.48 0.99998947 1.0000000 1.0000000 1.0000000
150 4.49 0.99998581 1.0000000 1.0000000 1.0000000
151 4.50 0.99998104 1.0000000 1.0000000 1.0000000
152 4.51 0.99997485 1.0000000 1.0000000 1.0000000
153 4.52 0.99996688 1.0000000 1.0000000 1.0000000
154 4.53 0.99995670 1.0000000 1.0000000 1.0000000
155 4.54 0.99994378 1.0000000 1.0000000 1.0000000
156 4.55 0.99992751 1.0000000 1.0000000 1.0000000
157 4.56 0.99990713 1.0000000 1.0000000 1.0000000
158 4.57 0.99988180 1.0000000 1.0000000 1.0000000
159 4.58 0.99985049 1.0000000 1.0000000 1.0000000
160 4.59 0.99981205 1.0000000 1.0000000 1.0000000
161 4.60 0.99976512 1.0000000 1.0000000 1.0000000
162 4.61 0.99970816 1.0000000 1.0000000 1.0000000
163 4.62 0.99963943 1.0000000 1.0000000 1.0000000
164 4.63 0.99955696 1.0000000 1.0000000 1.0000000
165 4.64 0.99945853 1.0000000 1.0000000 1.0000000
166 4.65 0.99934167 1.0000000 1.0000000 1.0000000
167 4.66 0.99920364 1.0000000 1.0000000 1.0000000
168 4.67 0.99904143 1.0000000 1.0000000 1.0000000
169 4.68 0.99885172 1.0000000 1.0000000 1.0000000
170 4.69 0.99863091 1.0000000 1.0000000 1.0000000
171 4.70 0.99837507 1.0000000 1.0000000 1.0000000
172 4.71 0.99808000 1.0000000 1.0000000 1.0000000
173 4.72 0.99774115 1.0000000 1.0000000 1.0000000
174 4.73 0.99735370 1.0000000 1.0000000 1.0000000
175 4.74 0.99691252 1.0000000 1.0000000 1.0000000
176 4.75 0.99641220 1.0000000 1.0000000 1.0000000
177 4.76 0.99584703 1.0000000 1.0000000 1.0000000
178 4.77 0.99521107 1.0000000 1.0000000 1.0000000
179 4.78 0.99449813 1.0000000 1.0000000 1.0000000
180 4.79 0.99370180 1.0000000 1.0000000 1.0000000
181 4.80 0.99281549 1.0000000 1.0000000 1.0000000
182 4.81 0.99183244 1.0000000 1.0000000 1.0000000
183 4.82 0.99074577 1.0000000 1.0000000 1.0000000
184 4.83 0.98954850 1.0000000 1.0000000 1.0000000
185 4.84 0.98823358 1.0000000 1.0000000 1.0000000
186 4.85 0.98679396 1.0000000 1.0000000 1.0000000
187 4.86 0.98522258 1.0000000 1.0000000 1.0000000
188 4.87 0.98351245 1.0000000 1.0000000 1.0000000
189 4.88 0.98165668 1.0000000 1.0000000 1.0000000
190 4.89 0.97964851 1.0000000 1.0000000 1.0000000
191 4.90 0.97748136 1.0000000 1.0000000 1.0000000
192 4.91 0.97514886 1.0000000 1.0000000 1.0000000
193 4.92 0.97264489 1.0000000 1.0000000 1.0000000
194 4.93 0.96996365 1.0000000 1.0000000 1.0000000
195 4.94 0.96709962 1.0000000 1.0000000 1.0000000
196 4.95 0.96404768 0.9999999 1.0000000 1.0000000
197 4.96 0.96080308 0.9999999 1.0000000 1.0000000
198 4.97 0.95736149 0.9999999 1.0000000 1.0000000
199 4.98 0.95371903 0.9999998 1.0000000 1.0000000
200 4.99 0.94987228 0.9999997 1.0000000 1.0000000
201 5.00 0.94581832 0.9999995 1.0000000 1.0000000
202 5.01 0.94155472 0.9999993 1.0000000 1.0000000
203 5.02 0.93707956 0.9999990 1.0000000 1.0000000
204 5.03 0.93239147 0.9999986 1.0000000 1.0000000
205 5.04 0.92748959 0.9999980 1.0000000 1.0000000
206 5.05 0.92237361 0.9999972 1.0000000 1.0000000
207 5.06 0.91704375 0.9999961 1.0000000 1.0000000
208 5.07 0.91150076 0.9999946 1.0000000 1.0000000
209 5.08 0.90574593 0.9999926 1.0000000 1.0000000
210 5.09 0.89978106 0.9999899 1.0000000 1.0000000
211 5.10 0.89360846 0.9999864 1.0000000 1.0000000
212 5.11 0.88723094 0.9999818 1.0000000 1.0000000
213 5.12 0.88065179 0.9999758 1.0000000 1.0000000
214 5.13 0.87387475 0.9999681 1.0000000 1.0000000
215 5.14 0.86690401 0.9999582 1.0000000 1.0000000
216 5.15 0.85974418 0.9999457 1.0000000 1.0000000
217 5.16 0.85240027 0.9999299 1.0000000 1.0000000
218 5.17 0.84487766 0.9999102 1.0000000 1.0000000
219 5.18 0.83718208 0.9998856 1.0000000 1.0000000
220 5.19 0.82931957 0.9998551 1.0000000 1.0000000
221 5.20 0.82129648 0.9998178 1.0000000 1.0000000
222 5.21 0.81311943 0.9997721 0.9999999 1.0000000
223 5.22 0.80479527 0.9997166 0.9999999 1.0000000
224 5.23 0.79633106 0.9996496 0.9999999 1.0000000
225 5.24 0.78773407 0.9995691 0.9999998 1.0000000
226 5.25 0.77901171 0.9994730 0.9999997 1.0000000
227 5.26 0.77017153 0.9993588 0.9999996 1.0000000
228 5.27 0.76122118 0.9992238 0.9999994 1.0000000
229 5.28 0.75216839 0.9990651 0.9999991 1.0000000
230 5.29 0.74302096 0.9988793 0.9999987 1.0000000
231 5.30 0.73378672 0.9986629 0.9999982 1.0000000
232 5.31 0.72447350 0.9984121 0.9999975 1.0000000
233 5.32 0.71508913 0.9981227 0.9999965 1.0000000
234 5.33 0.70564138 0.9977900 0.9999951 1.0000000
235 5.34 0.69613801 0.9974095 0.9999933 1.0000000
236 5.35 0.68658667 0.9969760 0.9999909 1.0000000
237 5.36 0.67699495 0.9964840 0.9999876 1.0000000
238 5.37 0.66737030 0.9959280 0.9999834 1.0000000
239 5.38 0.65772009 0.9953021 0.9999779 1.0000000
240 5.39 0.64805152 0.9946000 0.9999708 1.0000000
241 5.40 0.63837166 0.9938154 0.9999618 1.0000000
242 5.41 0.62868742 0.9929417 0.9999502 1.0000000
243 5.42 0.61900553 0.9919723 0.9999356 1.0000000
244 5.43 0.60933255 0.9909001 0.9999172 1.0000000
245 5.44 0.59967485 0.9897183 0.9998943 1.0000000
246 5.45 0.59003860 0.9884198 0.9998659 1.0000000
247 5.46 0.58042977 0.9869976 0.9998309 1.0000000
248 5.47 0.57085413 0.9854446 0.9997881 1.0000000
249 5.48 0.56131722 0.9837538 0.9997361 1.0000000
250 5.49 0.55182439 0.9819182 0.9996730 1.0000000
251 5.50 0.54238076 0.9799312 0.9995972 1.0000000
252 5.51 0.53299123 0.9777861 0.9995065 1.0000000
253 5.52 0.52366047 0.9754765 0.9993986 1.0000000
254 5.53 0.51439295 0.9729962 0.9992708 1.0000000
255 5.54 0.50519291 0.9703394 0.9991203 1.0000000
256 5.55 0.49606437 0.9675006 0.9989438 1.0000000
257 5.56 0.48701113 0.9644745 0.9987379 1.0000000
258 5.57 0.47803677 0.9612564 0.9984989 1.0000000
259 5.58 0.46914467 0.9578419 0.9982227 1.0000000
260 5.59 0.46033800 0.9542270 0.9979048 1.0000000
261 5.60 0.45161971 0.9504083 0.9975407 1.0000000
262 5.61 0.44299255 0.9463827 0.9971252 1.0000000
263 5.62 0.43445908 0.9421478 0.9966531 1.0000000
264 5.63 0.42602166 0.9377015 0.9961189 1.0000000
265 5.64 0.41768246 0.9330425 0.9955167 1.0000000
266 5.65 0.40944347 0.9281697 0.9948404 1.0000000
267 5.66 0.40130649 0.9230829 0.9940838 1.0000000
268 5.67 0.39327315 0.9177822 0.9932402 1.0000000
269 5.68 0.38534492 0.9122681 0.9923031 1.0000000
270 5.69 0.37752309 0.9065421 0.9912656 1.0000000
271 5.70 0.36980881 0.9006057 0.9901208 1.0000000
272 5.71 0.36220306 0.8944612 0.9888616 1.0000000
273 5.72 0.35470669 0.8881114 0.9874809 1.0000000
274 5.73 0.34732039 0.8815595 0.9859718 1.0000000
275 5.74 0.34004473 0.8748092 0.9843273 1.0000000
276 5.75 0.33288014 0.8678646 0.9825402 1.0000000
277 5.76 0.32582693 0.8607302 0.9806039 1.0000000
278 5.77 0.31888530 0.8534110 0.9785117 1.0000000
279 5.78 0.31205532 0.8459124 0.9762570 1.0000000
280 5.79 0.30533695 0.8382401 0.9738337 1.0000000
281 5.80 0.29873008 0.8303999 0.9712358 1.0000000
282 5.81 0.29223445 0.8223983 0.9684576 1.0000000
283 5.82 0.28584976 0.8142418 0.9654939 1.0000000
284 5.83 0.27957559 0.8059373 0.9623397 1.0000000
285 5.84 0.27341145 0.7974917 0.9589904 1.0000000
286 5.85 0.26735676 0.7889124 0.9554420 1.0000000
287 5.86 0.26141089 0.7802067 0.9516909 1.0000000
288 5.87 0.25557311 0.7713821 0.9477338 1.0000000
289 5.88 0.24984265 0.7624463 0.9435682 1.0000000
290 5.89 0.24421868 0.7534070 0.9391919 1.0000000
291 5.90 0.23870030 0.7442720 0.9346032 1.0000000
292 5.91 0.23328656 0.7350491 0.9298010 1.0000000
293 5.92 0.22797647 0.7257462 0.9247848 1.0000000
294 5.93 0.22276900 0.7163710 0.9195546 1.0000000
295 5.94 0.21766306 0.7069315 0.9141108 1.0000000
296 5.95 0.21265754 0.6974352 0.9084546 1.0000000
297 5.96 0.20775128 0.6878900 0.9025873 1.0000000
298 5.97 0.20294311 0.6783033 0.8965112 1.0000000
299 5.98 0.19823181 0.6686827 0.8902288 1.0000000
300 5.99 0.19361615 0.6590355 0.8837432 1.0000000
301 6.00 0.18909487 0.6493691 0.8770579 1.0000000
302 6.01 0.18466669 0.6396903 0.8701769 1.0000000
303 6.02 0.18033032 0.6300063 0.8631046 1.0000000
304 6.03 0.17608444 0.6203237 0.8558459 1.0000000
305 6.04 0.17192774 0.6106491 0.8484059 1.0000000
306 6.05 0.16785888 0.6009890 0.8407902 1.0000000
307 6.06 0.16387652 0.5913494 0.8330047 1.0000000
308 6.07 0.15997930 0.5817365 0.8250557 1.0000000
309 6.08 0.15616588 0.5721560 0.8169495 1.0000000
310 6.09 0.15243489 0.5626135 0.8086931 0.9999999
311 6.10 0.14878499 0.5531144 0.8002933 0.9999999
312 6.11 0.14521481 0.5436638 0.7917572 0.9999998
313 6.12 0.14172299 0.5342666 0.7830924 0.9999998
314 6.13 0.13830819 0.5249276 0.7743062 0.9999997
315 6.14 0.13496904 0.5156512 0.7654062 0.9999995
316 6.15 0.13170422 0.5064417 0.7564002 0.9999993
317 6.16 0.12851237 0.4973032 0.7472959 0.9999989
318 6.17 0.12539218 0.4882395 0.7381012 0.9999985
319 6.18 0.12234232 0.4792542 0.7288238 0.9999978
320 6.19 0.11936147 0.4703507 0.7194716 0.9999970
321 6.20 0.11644834 0.4615322 0.7100524 0.9999958
322 6.21 0.11360163 0.4528017 0.7005740 0.9999942
323 6.22 0.11082008 0.4441620 0.6910441 0.9999921
324 6.23 0.10810239 0.4356156 0.6814703 0.9999892
325 6.24 0.10544734 0.4271650 0.6718600 0.9999855
326 6.25 0.10285367 0.4188123 0.6622208 0.9999807
327 6.26 0.10032016 0.4105595 0.6525599 0.9999744
328 6.27 0.09784560 0.4024086 0.6428845 0.9999663
329 6.28 0.09542878 0.3943610 0.6332014 0.9999559
330 6.29 0.09306854 0.3864184 0.6235176 0.9999428
331 6.30 0.09076370 0.3785820 0.6138397 0.9999263
332 6.31 0.08851312 0.3708530 0.6041741 0.9999056
333 6.32 0.08631565 0.3632324 0.5945271 0.9998799
334 6.33 0.08417018 0.3557211 0.5849047 0.9998481
335 6.34 0.08207562 0.3483198 0.5753129 0.9998092
336 6.35 0.08003087 0.3410290 0.5657573 0.9997616
337 6.36 0.07803487 0.3338493 0.5562433 0.9997039
338 6.37 0.07608656 0.3267809 0.5467761 0.9996343
339 6.38 0.07418492 0.3198241 0.5373607 0.9995508
340 6.39 0.07232892 0.3129789 0.5280020 0.9994512
341 6.40 0.07051757 0.3062453 0.5187045 0.9993330
342 6.41 0.06874988 0.2996233 0.5094725 0.9991934
343 6.42 0.06702489 0.2931125 0.5003101 0.9990295
344 6.43 0.06534165 0.2867128 0.4912213 0.9988378
345 6.44 0.06369922 0.2804236 0.4822097 0.9986147
346 6.45 0.06209668 0.2742445 0.4732789 0.9983564
347 6.46 0.06053315 0.2681749 0.4644320 0.9980585
348 6.47 0.05900774 0.2622143 0.4556722 0.9977165
349 6.48 0.05751958 0.2563618 0.4470022 0.9973256
350 6.49 0.05606782 0.2506168 0.4384248 0.9968805
351 6.50 0.05465163 0.2449783 0.4299423 0.9963760Así, por ejemplo, la probabilidad estimada de que un sismo de magnitud superior a 6 ocurra en San Juan ocurra en los próximos 5 años es de casi \(0.65\), aumenta a \(0.88\) para los próximos 10 años y es de casi \(1\) para los próximos 20. En efecto, los únicos sismos registrados en los datos de tal magnitud son los de noviembre de 2016 y enero de 2021. Sismos de magnitud superior a 4.5 en San Juan son casi seguros para los próximos años, en todos los casos.
Código
knitr::kable(sismos_SJ %>%
filter(Magnitud > 6), caption = "Sismos de magnitud superior a 6 en San Juan, desde 2013 a 2021")| Fecha | Hora | Latitud | Longitud | Provincia | Percibido | Magnitud | Profundidad | FechaHora | Año |
|---|---|---|---|---|---|---|---|---|---|
| 2016-11-20 | 17:57:44 | -31.630 | -68.663 | San Juan | TRUE | 6.3 | 107 | 2016-11-20 17:57:44 | 2016 |
| 2021-01-18 | 23:46:20 | -31.854 | -68.963 | San Juan | TRUE | 6.4 | 8 | 2021-01-18 23:46:20 | 2021 |
5 Acerca de los datos
A continuación, se detallan aspectos de los datasets que conformaron el dataset reducido para el desarrollo del ejemplo, a la vez que se incluyen las fuentes de los datos y el código utilizado para pre-procesarlo con la sintaxis de tidyverse. De esta forma, puede fácilmente replicarse y/o adaptarse si así se lo desea.
También se incluye el enlace de descarga al dataset reducido, sismos-arg.csv, con el que se desarrolló el ejemplo.
5.1 Sobre el dataset de sismos
Los datos de este ejemplo corresponden a datos scrappeados del buscador de sismos del INPRES, cuyo buscador está disponible acá. En el dataset sismos_all.csv, disponible acá, en el que se registran eventos sísmicos en territorio argentino y diferentes variables asociadas a estos, en particular, si se trata de un sismo que se haya sentido o no.
Los sismos registrados en INPRES pueden estar revisados por un sismólogo y, posiblemente a razón de ello, algunas inconsistencias en la nomenclatura de las provincias es hallada. El siguiente toma el dataset de sismos_all.csv y:
- convierte en variable la percepción del sismo teniendo en cuenta que la intensidad solo se registra para sismos percibidos,
- limpia y unifica el formato de fecha de los eventos,
- extrae la información de profundidad, eliminando las unidades que forman parte del dato,
- filtra los datos para obtener los datos del territorio continental argentino (excluyendo a Tierra del Fuego),
- y recodifica los nombres de las provincias para unificar su identificación.
Código
# importo datos de sismos
sismos <- vroom("./fuente/04_visualizacion/sismos_all.csv")
sismos$Fecha <- as.Date(sismos$Fecha, tryFormats = c("%d/%m/%Y"))
# info de casi 10 años
min(sismos$Fecha)
max(sismos$Fecha)
# dummy para los percibidos
sismos <- sismos %>%
mutate(
Percibido = !is.na(Intensidad),
Magnitud = Magn.,
Profundidad = Profund.
)
# quito las unidades km. en intensidad
sismos$Profundidad <- as.numeric(str_sub(sismos$Profundidad, 1,-4))
# sismos en Argentina Continental (saco TdF porque se corre mucho el mapa)
sismos_arg <- sismos %>%
filter(
Provincia != "OCEANO PACIFICO",
Provincia != "CHILE",
Provincia != "MAULE - CHILE",
Provincia != "REPÚBLICA DE CHILE",
Provincia != "REPUBLICA DE CHILE",
Provincia != "SUR DE CHILE",
Provincia != "TFAIAS",
Provincia != "NORTH ISLAND, NEW ZEALAND",
Provincia != "REGION NOT FOUND.",
Provincia != "PERU",
Provincia != "PENINSULA ANTARTICA",
Provincia != "PASAJE DE DRAKE",
Provincia != "NORTHWEST OF KURIL ISLANDS",
Provincia != "NTARTIDA",
Provincia != "MAR DE SCOTIA",
Provincia != "MAR ARGENTINO",
Provincia != "ISLAS SANDWICH DEL SUR",
Provincia != "ISLAS SANDWICH",
Provincia != "ISLAS SHETLAND",
Provincia != "ISLAS SHETLAND DEL SUR",
Provincia != "ESTRECHO DE DRAKE",
Provincia != "FILIPINAS",
Provincia != "I.SANDWICH DEL SUR",
Provincia != "SECTOR ANTARTICO",
Provincia != "IS. SHETLAND DEL SUR",
Provincia != "ISLAS GEORGIA DEL SUR",
Provincia != "ISLAS GEORGIAS DEL SUR",
Provincia != "ISLAS GEORGIAS y SANDWICH DEL SUR",
Provincia != "ISLAS GIORGIA Y SANDWICH DEL SUR",
Provincia != "ISLAS ORCADAS",
Provincia != "ISLAS ORCADAS DEL SUR",
Provincia != "ATLANTICO SUR",
Provincia != "BOLIVIA",
Provincia != "PARAGUAY",
Provincia != "LIM. ARG-CHILE",
Provincia != "LIM.ARGENTINA-CHILE",
Provincia != "LIM.CHILE-ARGENTINA",
Provincia != "LIMITE ARGENTINA-CHILE",
Provincia != "LIMITE ARGENTINA CHILE",
Provincia != "TIERRA DEL FUEGO",
Provincia != "ISLAS GEORGIAS Y SANDWICH DEL SUR"
)
table(sismos_arg$Provincia)
# recodifico provincias
sismos_arg <- sismos_arg %>%
mutate(
Provincia = recode(
Provincia,
"BUENOS AIRES" = "Buenos Aires",
"CATAMARCA" = "Catamarca",
"CATAMARCA}" = "Catamarca",
"CATAMARCA LIM. TUCUMAN" = "Catamarca",
"CATAMARCA(LIM.CON TUCUMAN)" = "Catamarca",
"CATAMARCA(LIM.CON TUCUMÁN)" = "Catamarca",
"CHACO" = "Chaco",
"CHUBUT" = "Chubut",
"CORDOBA" = "Córdoba",
"CORDOBA(LIM.CON SAN LUIS)" = "Córdoba",
"CORDOBA LIM. SAN LUIS" = "Córdoba",
"CORRIENTES" = "Corrientes",
"ENTRE RIOS" = "Entre Ríos",
"FORMOSA" = "Formosa",
"jujuy" = "Jujuy",
"JUJUY" = "Jujuy",
"JUJUY(LIM.CON SALTA)" = "Jujuy",
"LA PAMPA" = "La Pampa",
"LA RIOJA" = "La Rioja",
"LA RIOJA LIM. SAN JUAN" = "La Rioja",
"LA RIOJA(LIM.CON CATAMARCA)" = "La Rioja",
"LA RIOJA(LIM.CON SAN JUAN)" = "La Rioja",
"LIM. CATAMARCA-TUCUMAN" = "Catamarca",
"LIM. LA RIOJA - SAN JUAN" = "La Rioja",
"LIMITE CATAMARCA - TUCUMAN" = "Catamarca",
"LIMITE CORDOBA - SAN LUIS" = "Córdoba",
"LIMITE JUJUY-SALTA" = "Jujuy",
"LIMITE JUJUY - SALTA" = "Jujuy",
"LIMITE LA RIOJA-CATAMARCA" = "La Rioja",
"LIMITE MENDOZA-SAN JUAN" = "Mendoza",
"LIMITE MENDOZA - SAN JUAN" = "Mendoza",
"LIMITE SALTA-JUJUY" = "Salta",
"LIMITE SALTA - CATAMARCA" = "Salta",
"LIMITE SALTA - JUJUY" = "Salta",
"LIMITE SAN JUAN-MENDOZA" = "San Juan",
"LIMITE SAN JUAN - LA RIOJA" = "San Juan",
"LIMITE SAN JUAN - MENDOZA" = "San Juan",
"LIMITE SAN JUAN SAN LUIS" = "San Juan",
"LIMITE SAN LUIS-SAN JUAN" = "San Luis",
"LIMITE TUCUMAN - CATAMARCA" = "Tucumán",
"MENDOZa" = "Mendoza",
"MENDOZA" = "Mendoza",
"MENDOZA(LIM.CON SAN JUAN)" = "Mendoza",
"MENDOZA}" = "Mendoza",
"MMENDOZA" = "Mendoza",
"NEUQUEN" = "Neuquén",
"RIO NEGRO" = "Río Negro",
"SALTA" = "Salta",
"SALTA (limite con Jujuy)" = "Salta",
"SALTA LIM. JUJUY" = "Salta",
"SALTA(LIM.CON JUJUY)" = "Salta",
"SALTA(LIM.CON TUCUMAN)" = "Salta",
"SAN JUAN LIM. MENDOZA" = "San Juan",
"SAN JUAN (LIM.CON SAN LUIS)" = "San Juan",
"SAN JUAN LIM. LA RIOJA" = "San Juan",
"SAN JUAN LIM.CON SAN LUIS" = "San Juan",
"SAN JUAN LIM.CON MENDOZA" = "San Juan",
"SAN JUAN LIM. MENDOZA" = "San Juan",
"SAN JUAN(LIM.ARG-CHI)" = "San Juan",
"SAN JUAN(LIM.CON LA RIOJA)" = "San Juan",
"LIMITE SAN JUAN - SAN LUIS" = "San Juan",
"LIMITE SAN JUAN LA RIOJA" = "San Juan",
"SAN JUAN(LIM.CON MENDOZA)" = "San Juan",
"SAN JUA" = "San Juan",
"SAn JUAN" = "San Juan",
"SAN JUAN" = "San Juan",
"LIMITE SAN JUAN MENDOZA" = "San Juan",
"SAN JUAN(LIM.CON SAN LUIS)" = "San Juan",
"SAN LUIS" = "San Luis",
"SANTA CRUZ" = "Santa Cruz",
"SANTIAGO DEL ESTERO" = "Santiago del Estero",
"SGO.DEL ESTERO LIM.CON CATAMARCA" = "Santiago del Estero",
"TUCUMAN" = "Tucumán",
"TUCUMÁN" = "Tucumán",
"TUCUMAN LIM. SALTA" = "Tucumán",
"TUCUMAN(LIM.CON CATAMARCA)" = "Tucumán"
)
)
sort(-table(sismos_arg$Provincia))
# Elijo variables
sismos_arg <- sismos_arg %>%
select(c(2, 3, 4, 5, 9, 10, 11, 12))5.1.1 Dataset pre-procesado: sismos-arg.csv
Para limpiar y reducir los datos al estudio de interés, se crea un dataset reducido, sismos-arg.csv, disponible acá, a partir de los datos scrappeados del buscador del INPRES. Se incluyen datos de sismos de diferente intensidad y magnitud en la región continental del país (excluyendo Tierra del Fuego) desde el 7 de enero de 2012, hasta el 18 de mayo de 2022, y se dispone de las siguientes variables.
Fecha: fecha del evento sísmico, en el formato año-mes-día.Hora: hora del evento sísmico, en el formato hora-minuto-segundo.Latitud: latitud del evento registrado.Longitud: longitud del evento registrado.Provincia: nombre de la provincia del evento registrado.Percibido:TRUEsi la magnitud fue percibida,FALSEsi no lo fue.Magnitud: magnitud del sismo en la escala de Richter.Profundidad: profundidad (km) registrada del sismo.
6 Librerías
Las librerías usadas para el desarrollo de este ejemplo, así como la información de la sesión de R, se muestran en el código que sigue.
Código
require(ggplot2)
require(tidyverse)
require(ggfortify)
require(plotly)
require(kableExtra)
require(vroom)
require(leaflet)
library(lubridate)Código
sessionInfo()R version 4.2.3 (2023-03-15 ucrt)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 10 x64 (build 22621)
Matrix products: default
locale:
[1] LC_COLLATE=Spanish_Argentina.utf8 LC_CTYPE=Spanish_Argentina.utf8
[3] LC_MONETARY=Spanish_Argentina.utf8 LC_NUMERIC=C
[5] LC_TIME=Spanish_Argentina.utf8
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] leaflet_2.1.2 vroom_1.6.1 kableExtra_1.3.4 plotly_4.10.1
[5] ggfortify_0.4.16 lubridate_1.9.2 forcats_1.0.0 stringr_1.5.0
[9] dplyr_1.1.2 purrr_1.0.1 readr_2.1.4 tidyr_1.3.0
[13] tibble_3.2.1 tidyverse_2.0.0 ggplot2_3.4.2
loaded via a namespace (and not attached):
[1] svglite_2.1.1 leaflet.providers_1.9.0 digest_0.6.31
[4] utf8_1.2.3 R6_2.5.1 evaluate_0.20
[7] httr_1.4.5 highr_0.10 pillar_1.9.0
[10] rlang_1.1.0 lazyeval_0.2.2 rstudioapi_0.14
[13] data.table_1.14.8 rmarkdown_2.21 labeling_0.4.2
[16] webshot_0.5.4 htmlwidgets_1.6.2 bit_4.0.5
[19] munsell_0.5.0 compiler_4.2.3 xfun_0.39
[22] pkgconfig_2.0.3 systemfonts_1.0.4 htmltools_0.5.5
[25] tidyselect_1.2.0 gridExtra_2.3 fansi_1.0.4
[28] viridisLite_0.4.2 crayon_1.5.2 tzdb_0.3.0
[31] withr_2.5.0 grid_4.2.3 jsonlite_1.8.4
[34] gtable_0.3.3 lifecycle_1.0.3 magrittr_2.0.3
[37] scales_1.2.1 cli_3.6.1 stringi_1.7.12
[40] farver_2.1.1 xml2_1.3.3 ellipsis_0.3.2
[43] generics_0.1.3 vctrs_0.6.1 tools_4.2.3
[46] bit64_4.0.5 glue_1.6.2 hms_1.1.3
[49] crosstalk_1.2.0 parallel_4.2.3 fastmap_1.1.1
[52] yaml_2.3.7 timechange_0.2.0 colorspace_2.1-0
[55] rvest_1.0.3 knitr_1.42